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Summary 

The frequency-phase dispersion and Pierce on-axis interac- 
tion impedance of a ferruled coupled-cavity traveling-wave- 
tube (TWT) slow-wave circuit were calculated using the 
three-dimensional simulation code Micro-SOS. The utiliza- 
tion of the code to reduce costly and time-consuming experi- 
mental cold tests is demonstrated by the accuracy achieved in 
calculating these parameters. A generalized input file was 
developed so that ferruled coupled-cavity TWT slow-wave 
circuits of arbitrary dimensions could be easily modeled. The 
practicality of the generalized input file was tested by apply- 
ing it to the ferruled coupled-cavity slow-wave circuit of the 
Hughes Aircraft Company model 961 HA TWT and by com- 
paring the results with experimental results. 

Introduction 

Micro-SOS is a three-dimensional electromagnetic circuit 
analysis code developed by Mission Research Corporation 
(MRC) in collaboration with Varian Associates as part of a 
project sponsored by the Department of Defense (DOD) and 
the National Aeronautics and Space Administration (NASA) 
to develop computer software for microwave tube design 
(refs. 1 to 3). Micro-SOS is a version of a larger code, SOS 
(self-optimized sector), which includes the simulation of 
charged-particle trajectories. A frequency-domain solver was 
added to SOS and the particle-handling ability was deleted to 
form Micro-SOS. The finite-difference code is written in 
FORTRAN 77 and can be operated in either a time-domain 
algorithm in which Maxwell’s time-dependent equations are 
solved or in a frequency-domain algorithm in which the reso- 
nant mode frequencies and field patterns are obtained. The 
results described in this report were obtained entirely from 
simulations using the frequency-domain algorithm of Micro- 
SOS. 

The objectives of the code development were to provide an 
electromagnetic circuit modeling/analysis capability that was 
better, faster, and cheaper than that previously available and 
to provide a means of obtaining more information about the 
radiofrequency (RF) fields in a circuit than is normally pos- 
sible through experimental hardware “cold testing.” Hard- 
ware cold tests refer to evaluating circuit models (actual or 
scaled) on the RF test bench with no electron beam present. 


The cold tests are used to determine the RF phase shift, inter- 
action impedance, and attenuation for a cavity. These data 
can then be used as input for a coupled-cavity TWT RF- 
beam-interaction computer model (ref. 4) which simulates 
the output RF power and beam characteristics. 

Using Micro-SOS instead of experimental cold testing to 
obtain these data has several advantages: significant conser- 
vation of labor, time, and cost; efficient exploration of novel 
circuit structure designs; determination of higher frequency 
mode information which avoids expensive experimental tests 
to determine potential RF instabilities in a given circuit 
(ref. 5). 

The practicality and efficiency of the frequency-domain 
algorithm of Micro-SOS were demonstrated by simulating 
the frequency-phase dispersion and Pierce on-axis interaction 
impedance characteristics for a ferruled coupled-cavity slow- 
wave circuit. The Micro-SOS results obtained for both the 
dispersion and impedance characteristics are compared in this 
report to experimental cold-test results. The facility of using 
the generalized Micro-SOS input file for ferruled coupled- 
cavity TWT circuits is also shown. By simply replacing the 
existing dimensions with the dimensions of the circuit of in- 
terest, dispersion curves and the Pierce on-axis interaction 
impedance characteristics can be readily determined for any 
ferruled coupled-cavity TWT circuit. 

Symbols 

a n ,b n ,E n on-axis axial electric field harmonic magnitude, V/m 
E ztot total on-axis axial electric field, V/m 
/ frequency, GHz 

g ferrule gap, m 

h cavity gap, m 

K n Pierce on-axis interaction impedance, ft 
L cavity length, m 

N number of cavities in simulation 

n space harmonic order 
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P period of E^, m 

Prf RF power flow, J/s 

r\ ferrule inner radius, m 

r 2 ferrule outer radius, m 

/*3 slot inner radius, m 

r 4 slot outer radius, m 

r$ cavity radius, m 

t cavity wall thickness, m 

v g group velocity, m/s 

W electromagnetic energy per unit length, J/m 

Wf total electromagnetic energy, J 

z axial coordinate, m 

P axial phase constant, rad/m 

0 slot arc, deg 

< p phase shift per cavity for fundamental harmonic, 

rad 

ft) angular frequency, rad 

Background 

Micro-SOS was used to model the ferruled coupled-cavity 
slow- wave circuit of the Hughes Aircraft Company 961 HA 
TWT currently under development (ref. 6). The 961 HA, a 
59- to 64-GHz TWT for intersatellite communications, was 
designed to have an overall efficiency exceeding 40 percent 
over the 5-GHz bandwidth for dual-saturated RF output 
power levels of 30 and 75 W. The design is a revision of the 
earlier model, the 96 1H (ref. 7), which had a slot-mode oscil- 
lation problem caused by the interaction of the electron beam 
with the slot mode at the upper cutoff frequency. When the 
dimensions of the coupling slot in the 961 HA were adjusted 
to shift the slot mode down in frequency, this problem was 
resolved (ref. 7). 

Two views of the generalized ferruled coupled-cavity cir- 
cuit are shown schematically in figures 1(a) and (b). The cou- 
pling slots are rotated 180° at alternating cavity partition 
walls. The ferrules, which are the hollow posts surrounding 


the beam hole, concentrate the RF electric field in the beam 
region to increase the interaction impedance of the cavity. 
The corresponding dimensions of a single cavity of the 
961HATWT are listed in table I. The circuit implements 169 
cavities over 3 sections with a mild two-stage velocity taper 
in the output section (ref. 6). 


TABLE I.— DIMENSIONS OF 961 HA 
FERRULED COUPLED-CAVITY 
TRAVELING- WAVE-TUBE 
CIRCUIT 

[See fig. L] 


Ferrule radius, 10* 3 m 

Inner, r x 0.2794 

Outer, r 2 0.4318 

Slot radius, 10" 3 m 

Inner, r 3 0.6858 

Outer r 4 1.1049 

Cavity radius, r 3 , 1 0 -3 m 1.1 049 

Gap, 10" 3 m 

Ferrule, g 0.2921 

Cavity,/! 0.7518 

Wall thickness, f, 1 0~ 3 m 0.2 1 59 

Cavity length, L, 10“ 3 m 0.9677 

Slot arc, 0, deg 140 


Analysis 

Dispersion Simulation 

The frequency-phase dispersion characteristics of the fer- 
ruled coupled-cavity slow-wave circuit were obtained by 
modeling one- to four-cavity configurations of the structure. 
The frequency algorithm of Micro-SOS calculated the reso- 
nant frequencies for different values of phase shift per cavity 
to form an accurate dispersion curve. The different eigen- 
modes are forced by axially truncating the different cavity 
configurations with either an electric wall (simulated by using 
the STRUCTURE command) or a magnetic wall (simulated 
by using the SYMMETRY MIRROR command). An electric 
wall is a boundary that simulates a perfect conductor with the 
electric field perpendicular to the wall. At a magnetic wall 
boundary condition, the tangential magnetic field is zero and 
the magnetic field is perpendicular to the wall. 

Truncating the longitudinal ends of the circuit with these 
boundary conditions corresponds to simulating standing 
waves with an integral number of half-wavelengths (phase 
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(b) Top view. 

Figurel . — Generalized ferruled coupled-cavity, traveling wave 
tube (TWT) circuit. 


shifts of Ji) in the length of the circuit. The cavity configura- 
tions and boundary conditions required for obtaining several 
resonant frequencies for the first forward RF space harmonic 
of the cavity mode and the first backward RF space harmonic 
of the slot mode using Micro-SOS are listed in table II 
(ref. 3). For each combination of cavity configuration and 
boundary conditions, the phase shifts per cavity for the 
resulting resonance frequencies are listed in the order of 
increasing frequency, with the slot-mode phase shifts 
enclosed in parentheses. 

The total cavity-mode axial on-axis electric field E ztot is 
shown in figure 2 for each of the phase shifts per cavity listed 
in table II. The plots are extended over 12 cavities in order to 
show the full periodic nature of each case. The locations of the 
boundaries in the Micro-SOS simulations are shown by the 


TABLE II.— BOUNDARY CONDITIONS FOR 
RESONANCE AT VARIOUS PHASE 
SHIFTS PER CAVITY 


Number of 
cavities 

Boundary 

conditions 3 

Phase shift per cavity, b 

PL 

rad 

1 

MM 

71 

(2 n, n) 

1 

ME 

3t*2 

(3a/2) 

1 

EE 

2n 

2 

ME 

5t*4, 7 7*4 
(7 7*4, 5 ti/4) 

2 

EE 

3t*2, 2tt 

3 

ME 

Ijdb.lidl, 11 ti/6 
(IIji/6, 3n/2, 7 te/6) 

3 

EE 

4t*3, 5?*3, 2 7T 
(57*3, 4t*3) 

4 

EE 

5 ti/ 4, 3 n/2, 7 ti/4, lit 
(771/4, 3ti/2, 5ti/4) 


"Magnetic wall, M; electric wall, E. 
b Slot-mode phase shifts are in parentheses. 


vertical truncation lines in figure 2, with electric and magnetic 
walls indicated by solid and dashed lines, respectively. 

Because of the relatively small coupling slot size in this 
coupled-cavity circuit, the slot mode (which is technically a 
separate passband in the same mode as the cavity-mode pass- 
band) is at a higher frequency than the cavity mode (refs. 5 
and 8). The axial electric field distribution near the axis is 
similar to that of the cavity mode at the same value of phase 
shift per cavity, except at 2k where the axial electric field is 
zero at the axis. The slot mode is characterized by a strong 
radial component of the electric field in some or all of the 
coupling slots. 

A generalized input file was developed so that the resonant 
frequencies could be determined for any ferruled coupled- 
cavity TWT slow-wave circuit. As listed in the appendix, the 
file simulates four cavities with two electric walls. However, 
the file can be easily modified to simulate one, two, or three 
cavities with any of the three boundary condition combina- 
tions. The number of cavities modeled N is defined by the 
user in the CAVITY NUMBER section of the appendix; the 
boundary conditions can be adjusted by editing the SYMME- 
TRY section. 

In the generalized input file listed in the appendix, the grid 
is defined with the X1GRID, X2GRID, and X3GRID com- 
mands which refer to the radial, azimuthal, and axial direc- 
tions, respectively. To model general ferruled coupled-cavity 
circuits in addition to that of the model 961 HA TWT of this 
report, DEFINE commands are used frequently throughout 
the input file of the appendix. By exercising the DEFINE 
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Phase shift 
per cavity, 



Figure 2. — Total axial on-axis electric field over 1 2 cavities with in- 
dicated boundary conditions (magnetic wall, M; electric wall, E). 


command option, it is possible to specify the entire input file 
in terms of variable names. The CONVERT command used 
in the PARAMETER section of the appendix converts the 
dimensions of the cavity from inches to meters and simulta- 
neously assigns that dimension to a variable name. If the user 
wishes to enter the dimensions originally in meters, the 
CONVERT commands can be replaced with DEFINE com- 
mands (ref. 2). 

This generalized input file allows the user some freedom 
to experiment with various dimensions by changing 
the parameter dimensions in the PARAMETER section. It 
also allows experimentation with the grid spacings. Serving 
as a guide for the GRID SPACING section, figures 3(a) to (c) 
show views of one cavity of the ferruled coupled-cavity 
structure with the xl-, x2-, and x3-directions divided into 
regions. In this section, one can select a grid spacing for each 
of these regions. The file gives default values which have 
been proven to provide accurate results for the slow-wave 
circuit of the 961HA TWT. 


X1REG1 — ► 

X1REG2 — ► 

X1REG3 — ► 
X1REG4 — ► 



mmmz. 





(c) Axial direction x3. 

Figure 3. — Generalized ferruled coupled- cavity TWT 
circuit simulation segments. 
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Impedance Simulation 


with 


For each value of phase shift per cavity on the dispersion 
curve, the Pierce on-axis interaction impedance for the n th RF 
space harmonic (ref. 9) is defined as 


K 


n 


^Pn P RF 


( 1 ) 


P 

2 2 r 

a „=7 } E zto, COS {Pn Z ) dz {1) 

P 

~2 


th 

where E n is the on-axis electric field magnitude of the n 
harmonic, is the axial phase constant of the n th harmonic 
defined by 

(2) 

where <p is the phase shift per cavity in radians of the funda- 
mental harmonic, n is the space harmonic order, L is the 
length of one cavity, and P^f is the RF power flow defined 
by 

P RF = Wv g (?) 


and 

P 

2 r 

b n=J (8) 

_P 

where P is the period of E ztot . By defining z = 0 at a bound- 
ary, E ztot will always be either an even or an odd function 
with respect to z. If z = 0 at an electric wall, E ztot is an even 
function with b n = 0, and the magnitude of the first harmonic 
is obtained from equation (7): 


where is the group velocity, and W is the stored electro- 
magnetic energy per unit length, 



(4) 


with the total energy Wj calculated with Micro-SOS in the 
number of cavities simulated N. Thus, the Pierce on-axis in- 
teraction impedance can be expressed as 


, E%NL 
" " 2p 2 n W T v x 


(5) 


P 

2 

E \ = j\ ‘ E zto, cos M dz (9) 

o 

If z = 0 at a magnetic wall, E ztot is an odd function with 
a n - 0, and the magnitude of the first harmonic is obtained 
from equation (8): 


P 

E \ = jj E z'oMPl Z ) dz (10) 

0 


In the typical ferruled coupled-cavity TWT circuit, the 
beam is synchronous only with the first RF space harmonic. 
Thus, we need to calculate only the first (n = 1) harmonic 
term in equation (5). To determine its magnitude E\, a Fou- 
rier analysis is performed on the total on-axis axial electric 
field E ztoh which is tabulated at the center of each grid space 
along the axis in the z-direction by using the LINPRINT 
command. 

In the analysis, the total axial electric field on the axis can 
be represented in terms of a Fourier series, 


The group velocity of equation (5) is defined as 


v 


8 


do) 


(in 


and is extracted from the dispersion curve which is plotted 
with frequency/(in gigahertz given by Micro-SOS) versus 
p L (in radians). Using these parameters, y*, can be expressed 
as 


E z'o, =X a « COS(/3 ^ ) + fo « Sin (^ Z ) (6) 


V =2 7tL , 

* 


( 12 ) 
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Finally, to obtain Wj for equation (5), the BALANCE 
command is used to control the Micro-SOS energy calcula- 
tion. The OBSERVE command is then used to record the 
specified variables during simulation. By specifying EN- 
ERGY and TOTAL in the OBSERVE command and desig- 
nating the entire volume of the structure, the total electro- 
magnetic energy Wfis obtained. 


Results 

The ferruled coupled-cavity slow-wave circuit of the 
model 961 HA TWT was modeled in both the Cartesian and 
cylindrical coordinate systems. The results from these simu- 
lations were then compared with experimental results. 

Dispersion Simulations 

The Cartesian model was developed first. The Aug89SOS 
version of Micro-SOS (refs.l to 3) was used to generate a 
grid cell resolution of 39 by 39 in the transverse plane and 1 3 
per cavity in the longitudinal direction. Figures 4(a) and (b), 
which were generated by Micro-SOS, show the grid only in 
the conductor regions. The lines of the mesh were spaced as 
equally as possible and grid boundaries were matched with 
major geometrical features such as the beam hole and fer- 
rules. The resonant frequencies obtained by modeling one 
and two cavities are compared with experimental data for the 
cavity mode in figure 5 and in table III. 


TABLE III.— CAVITY-MODE RESONANT 
FREQUENCIES FOR 961 HA USING 
CARTESIAN COORDINATES 


Phase shift 
per cavity, 

PL, 

rad 

Experimental 

frequency, 

/ 

GHz 

Micro-SOS 

frequency, 

/. 

GHz 

Frequency 

difference 

v, 

percent 

71 



52.53 



5ti/4 

55.88 

54.67 

-2.17 

3ti/2 

61.38 

60.35 

-1.68 

ln/4 

69.87 

68.30 

-2.25 

2n 

75.47 

74.48 

-1.31 


In comparison with data provided by the manufacturer 
(H.C. Limburg, 1992, Hughes Aircraft Co., Torrance, CA, 
personal communication), the Micro-SOS results using the 
Cartesian model were consistently lower by an average of 
1.85 percent. Increasing the resolution of the grid resulted in 
only a slight reduction of the difference at the cost of in- 
creased computation time. It should be noted that in reference 
3 the calculated resonance frequencies were reported to be 
consistently lower than the experimental results by about 4 
percent for a similar ferruled coupled-cavity TWT 



(a) End view. 



x3 

(b) Top view. 


Figure 4. — Micro-SOS grid for 961 HA ferruled coupled-cavity 
TWT circuit over two cavities. 


To determine if more accurate results could be obtained, 
the ferruled coupled-cavity slow-wave circuit was also mod- 
eled in cylindrical coordinates with a grid resolution of 13 by 
18 by 13 per cavity in the radial, azimuthal, and axial direc- 
tions, respectively. Again, the grid lines were spaced as 
equally as possible while the grid boundaries were kept coin- 
cident with all surfaces. 
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Figure 5. — Experimental and Micro-SOS Cartesian coordinate 
simulation of dispersion for cavity mode of 961 HA TWT. 


Results for the dispersion curve simulated with one-, 
two-, and three-cavity configurations in the cylindrical coor- 
dinate system are shown in figure 6 and are presented in 
tables IV and V for the cavity and slot modes. Compared 
with the experimental data, the cavity-mode results agree 
within 0.37 percent for the lower frequency resonances of the 

Mode 

O Experimental Cavity 

♦ Micro-SOS Cavity 

O Experimental Slot 



Figure 6. — Experimental and Micro-SOS cylindrical coordinate 
simulations of dispersion for cavity and slot modes of 961 HA 
TWT. 


TABLE IV — CAVITY-MODE RESONANT 
FREQUENCIES FOR 961 HA TRAVELING 
WAVE TUBE USING CYLINDRICAL 
COORDINATES 


Phase shift 
per cavity, 

PL, 

rad 

Experimental 

frequency, 

L 

GHz 

Micro-SOS 

frequency, 

/. 

GHz 

Frequency 

difference 

v. 

percent 

K 

— 

53.69 

— 

ln/6 

54.44 

54.64 

+0.37 

5t/4 

55.88 

55.74 

-.25 

371/2 

61.38 

61.18 

-.33 

ln/4 

69.87 

68.59 

-1.83 

11 71/6 

72.03 

71.00 

-1.43 

2tt 

75.47 

73.68 

-2.37 


TABLE V — SLOT-MODE RESONANT 
FREQUENCIES FOR 961 HA TRAVELING 
WAVE TUBE USING CYLINDRICAL 
COORDINATES 


Phase shift 
per cavity, 
PL, 

rad 

Experimental 

frequency, 

l 

GHz 

Micro-SOS 

frequency, 

/, 

GHz 

Frequency 

difference 

v. 

percent 

TC 



112.83 

— 

ln/6 

112.71 

111.08 

-1.45 

5t/4 

110.91 

109.37 

-1.39 

3t/2 

102.98 

101.81 

-1.14 

In! 4 

93.25 

92.86 

-0.42 

11 71/6 

90.30 

90.16 

-0.16 

2n 

— 

87.72 

— 


dispersion curve, but for increasing frequencies, the differ- 
ences increase significantly to -2.37 percent for a phase shift 
per cavity of In. The Micro-SOS cylindrical results were sig- 
nificantly closer to the experimental results with an average 
absolute frequency difference of only 1.10 percent, compared 
to the Cartesian results with the average difference being 
1.85 percent. Increasing the grid resolution did not signifi- 
cantly reduce the difference. The Micro-SOS results for the 
slot mode agree very well with the experimental data with an 
average absolute frequency difference of 0.91 percent. 

Impedance Simulation 

The Pierce on-axis interaction impedance for the first for- 
ward space harmonic was calculated from equation (5) using 
the Micro-SOS cylindrical model dispersion simulation 
results. The impedance results are shown in figure 7 and 
table VI. As shown, the impedance values match very well, 
especially for the lower values of phase shift, with an average 
absolute difference of 3.09 percent when compared 
with experimental cold-test results. 
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TABLE VI — PIERCE ON-AXIS INTERACTION 
IMPEDANCE 


Phase shift 
per cavity, 

PL, 

rad 

Experimental 

impedance, 

*i. 

Micro-SOS 

impedance, 

Q 

Impedance 

difference 

AK, 

percent 

5jl/4 

23.85 

23.83 

-0.08 

4tc/3 

14.56 

14.62 

0.41 

3n/2 

7.09 

7.32 

3.24 

5jc/3 

4.26 

4.59 

7.75 

7jc/4 

3.53 

3.67 

3.97 



Figure 7. — Experimental and Micro-SOS cylindrical coordinate 
simulation of Pierce on -axis interaction impedance for 
961HATWT. 


Conclusions 

A generalized ferruled coupled-cavity traveling-wave-tube 
(TWT) circuit input file was developed for the three-dimen- 
sional electromagnetic circuit analysis code Micro-SOS. Ar- 
bitrary circuit geometries can be easily modeled by simply 
adjusting dimensions in the file. The use of the file for the 
accurate generation of the frequency-phase dispersion and 
Pierce on-axis interaction impedance characteristics was 
detailed. 

The generalized input file was used to simulate the disper- 
sion and impedance for a specific circuit, that of the 961 HA 
coupled-activity TWT The agreement with experiment was 
very good with an average absolute frequency difference of 
about 1 percent and an average absolute impedance differ- 
ence of about 3 percent. These results demonstrate the reli- 
ability of Micro-SOS simulations and their capability to re- 
duce expensive and time-consuming experimental procedures 
in the design process for coupled-cavity TWT's. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio 
October 4, 1993 
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APPENDIX— Micro-SOS Input File 


TITLE * GENERALIZED INPUT FILE FOR A FERRULED COUPLED-CAVITY CIRCUIT 
IN CYLINDRICAL COORDINATES * / 

SYSTEM CYLINDRICAL / 


0 ************************************************************************* / 

C PARAMETER SECTION PARAMETER SECTION PARAMETER SECTION / 

Q ************************************************************************* / 

COMMENT * CONVERT AND DEFINE DIMENSIONS, DIMENSIONS LISTED ARE FOR THE 

HUGHES 961HA CIRCUIT * / 


COMMENT * FERRULE INNER RADIUS * rl * / 

CONVERT RBEAMHOLE 0.0110 INCHES t 

COMMENT * FERRULE OUTER RADIUS * r2 * / 

CONVERT RFERRULE 0.0170 INCHES / 

COMMENT * SLOT INNER RADIUS «= r3 * / 

CONVERT RPATCH 0.0270 INCHES / 

COMMENT * SLOT OUTER RADIUS = r4 * / 

CONVERT RSLOT 0.0435 INCHES / 

COMMENT * CAVITY RADIUS - r5 * / 

CONVERT RCAVITY 0.0435 INCHES / 

COMMENT * FERRULE GAP “ g * / 

CONVERT FERRULESPACE 0.0115 INCHES / 

COMMENT * CAVITY GAP "= h * / 

CONVERT CAVITYSPACE 0.0296 INCHES t 

COMMENT * WALL THICKNESS * t * / 

CONVERT WALLTHICK 0.0085 INCHES / 

COMMENT * PERIOD LENGTH *= L * / 

CONVERT CAVITYLENGTH 0.0381 INCHES / 

COMMENT * SLOT ARC = theta * / 

DEFINE SLOTARC 140 / 

♦ 

DEFINE PI 3.14159265359 / 

0 A************************************************************************* / 

0 **************************** *********************************************** / 


0 *★★★*★★★★*★★**★★★★****★*★*★**★*★★*★•*****★★***★*★*★******★******★*★****★**★ / 
C CAVITY NUMBER SECTION. . .CAVITY NUMBER SECTION. . .CAVITY NUMBER SECTION. . . 

0 ************************************************************************** / 
COMMENT * DEFINE THE NUMBER OF CAVITIES TO BE MODELED, N * / 

DEFINE N 4 / 

0 a************************************************************************* / 

0 a************************************************************** *********** f 

0 a**************************************************************'*********** / 

C GRID SPACING SECTION GRID SPACING SECTION GRID SPACING SECTION / 

0 a********************** a***************************************** ********* / 

COMMENT * FOLLOW FIGURES 3a, b AND c FOR REGIONS * / 

DEFINE X1REG1 4 / 

DEFINE X1REG2 3 / 

DEFINE X1REG3 2 / 

DEFINE X1REG4 4 / 



DEFINE X2REG1 1 / 
DEFINE X2REG2 7 / 


DEFINE X3REG1 1 / 
DEFINE X3REG2 3 / 
DEFINE X3REG3 5 / 


DEFINE CAVSTART CAVITYLENGTH CAVITYSPACE SUBTRACT 2.0 DIVIDE / 

DEFINE LENGTH CAVITYLENGTH N MULTIPLY / 

DEFINE X1MAX RCAVITY / 

DEFINE XIMIN 0 / 

DEFINE X2MAX PI 2.0 MULTIPLY / 

DEFINE X2MIN 0 / 

DEFINE X3MAX LENGTH 2 DIVIDE / 

DEFINE X3MIN X3MAX -1 MULTIPLY / 

DEFINE X1TOTGRID X1REG1 X1REG2 ADD X1REG3 ADD X1REG4 ADD 1 ADD / 

DEFINE X2TOTGRID X2REG1 4 MULTIPLY X2REG2 ADD X2REG2 ADD 1 ADD / 

DEFINE X3PER X3REG1 X3REG2 ADD X3REG3 ADD X3REG2 ADD X3REG1 ADD / 

DEFINE X3TOTGRID X3PER N MULTIPLY 1 ADD / 

DEFINE XI BLOCK X1TOTGRID 1 ADD / 

DEFINE X2BLOCK X2TOTGRID 1 ADD / 

DEFINE X3BLOCK X3TOTGRID 1 ADD / 

DEFINE X1SPACE1 RCAVITY RPATCH MINUS X1REG1 DIVIDE / 

DEFINE X1SPACE2 RPATCH RFERRULE MINUS X1REG2 DIVIDE / 

DEFINE X1SPACE3 RFERRULE RBEAMHOLE MINUS X1REG3 DIVIDE / 

DEFINE X1SPACE4 RBEAMHOLE X1REG4 DIVIDE / 

DEFINE X1DIM1 X1SPACE1 X1REG1 MULTIPLY / 

DEFINE X1DIM2 X1SPACE2 X1REG2 MULTIPLY / 

DEFINE X1DIM3 X1SPACE3 X1REG3 MULTIPLY / 

DEFINE X1DIM4 X1SPACE4 X1REG4 MULTIPLY / 

DEFINE X2DIM2 SLOTARC PI MULTIPLY 180 DIVIDE / 

DEFINE X2DIM1 PI X2DIM2 MINUS 2 DIVIDE / 

DEFINE X2SPACE1 X2DIM1 X2REG1 DIVIDE / 

DEFINE X2SPACE2 X2DIM2 X2REG2 DIVIDE / 

DEFINE X3SPACE1 CAVSTART X3REG1 DIVIDE / 

DEFINE FERRIIEIGHT CAVITYLENGTH FERRULESPACE MINUS 2 DIVIDE CAVSTART MINUS / 
DEFINE X3SPACE2 FERRHEIGHT X3REG2 DIVIDE / 

DEFINE X3SPACE3 FERRULESPACE X3REG3 DIVIDE / 

DEFINE X3DIM1 CAVSTART / 

DEFINE X3DIM2 X3SPACE2 X3REG2 MULTIPLY / 

DEFINE X3DIM3 X3SPACE3 X3REG3 MULTIPLY / 


X1GRID 

FUNCTION 

X1TOTGRID 

1 

XIMIN 

X1REG4 

X1SPACE4 

X1DIM4 






X1REG3 

XISPACE3 

X1DIM3 






X1REG2 

X1SPACE2 

X1DIM2 






X1REG1 

X1SPACE1 

X1DIM1 / 

X2GRID 

FUNCTION 

X2TOTGRID 

1 

X2MIN 

X2REG1 

X2SPACE1 

X2DIM1 






X2REG2 

X2SPACE2 

X2DIM2 






X2REG1 

X2SPACE1 

X2DIM1 






X2REG1 

X2SPACE1 

X2DIM1 






X2REG2 

X2SPACE2 

X2DIM2 
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X2REG1 

X2SPACE1 

X2DIM1 

/ 

X3GRID FUNCTION X3TOTGRID 1 X3MIN X3REG1 

X 3 SPACE 1 

X3DIM1 


X3REG2 

X3SPACE2 

X3DIM2 


X3REG3 

X3SPACE3 

X3DIM3 


X3REG2 

X3SPACE2 

X3DIM2 


X3REG1 

X3SPACE1 

X3DIM1 


X3REG1 

X3SPACE1 

X3DIM1 


X3REG2 

X3SPACE2 

X3DIM2 


X3REG3 

X3SPACE3 

X3D1M3 


X3REG2 

X3SPACE2 

X3DIM2 


X3REG1 

X 3 SPACE 1 

X3DIM1 


X3REG1 

X3SPACE1 

X3D1M1 


X3REG2 

X3SPACE2 

X3DIM2 


X3REG3 

X3SPACE3 

X3DIM3 


X3REG2 

X3SPACE2 

X3DIM2 


X3REG1 

X3SPACE1 

X3DIM1 


X3REG1 

X3SPACE1 

X3DIM1 


X3REG2 

X3SPACE2 

X3DIM2 


X3REG3 

X3SPACE3 

X3DIM3 


X3REG2 

X3SPACE2 

X3DIM2 


X3REG1 

X3SPACE1 

X3DIM1 

/ 


Q ****************************************************************** ******** / 

C SYMMETRY SECTION SYMMETRY SECTION SYMMETRY SECTION / 

C ************************************************************************** / 

SYMMETRY AXIAL *1111 X2TOTGRID 1 X3TOTGRID / 

SYMMETRY PERIODIC +1 1 X1TOTGRID 11 1 X3TOTGRID 

ANTI-ALIGN 1 X1TOTGRID X2TOTGRID X2TOTGRI D 1 X3T0TGRID / 


COMMENT * SEE TABLE II. FOR BOUNDARY/CAVITY CONFIGURATIONS. 

2 SYMMETRY PLANES SIMULATE THE MM BOUNDARY CONDITION, 

COMMENTING OUT 1 SYMMETRY SIMULATES AN EM BOUNDARY CONDITION, 
COMMENTING OUT BOTH SYMMETRY PLANES SIMULATES AN EE BOUNDARY * / 

C SYMMETRY MIRROR ALIGN 1 X1T0TGRID 1 X2TOTGRID 11 / 

C SYMMETRY MIRROR ANTI-ALIGN 1 X1TOTGRID 1 X2TOTGRID X3TOTGRID X3TOTGRID / 

C ************************************************************************* / 

Q ************************************************************************* / 

STRUCTURE CONFORMAL SOLIDIFY +1 1 X1BLOCK 1 X2TOTGRID 0 X3BLOCK / 

STRUCTURE CONFORMAL CAVITY -1 1 X1TOTGRID 1 X2TOTGRID 1 X3TOTGRID / 

DEFINE IFERRST X1REG4 1 ADD / 

DEFINE IFERREN IFERRST X1REG3 ADD / 

STRUCTURE CONFORMAL FERRULE +1 IFERRST IFERREN 1 X2TOTGRID 1 X3TOTGRID / 

£» **************★***** + **★★****★*★**★**** + ***★***★★★★*♦*♦♦*★**■***♦*♦****** / 
C WALL SECTION . . .WALL SECTION. . .WALL SECTION. . .WALL SECTION. . . / 

c ************************************************************************ / 

COMMENT * FORM THE WALLS * / 

DEFINE LEND X3REG1 1 ADD / 

STRUCTURE CONFORMAL LWALL +1 IFERREN X1TOTGRID 1 X2TOTGRID 1 LEND / 

DEFINE RST X3TOTGRID X3REG1 MINUS / 

STRUCTURE CONFORMAL RWALL +1 IFERREN X1TOTGRID 1 X2TOTGRID RST X3TOTGRID / 


DEFINE N1 N 1 MINUS / 

COMMENT * FORM THE CENTER WALLS * / 
DO K 1 Ni 1 / 



DEFINE M1END LEND X3PER K MULTIPLY ADD / 

DEFINE MIST M1END X3REG1 MINUS X3REG1 MINUS / 

STRUCTURE CONFORMAL MWALLN +1 IFERREN X1TOTGRID 1 X2TOTGRID MIST M1END / 

ENDDO / 

C A********************************************************************* / 

Q ****★*******★***********★★**************★**★***•*•***★**★★★*****★★★***** j 

Q j 

C SLOT SECTION ... SLOT SECTION ... SLOT SECTION. . -SLOT SECTION... / 

C ********************************************************************** / 

COMMENT * FORM THE SLOTS * / 


DEFINE ISLOTST IFERREN X1REG2 PLUS / 

DEFINE SLRST X2REG1 1 ADD / 

DEFINE SLREND SLRST X2REG2 ADD / 

DEFINE SLLST SLREND X2REGI ADD X2REG1 ADD / 

DEFINE SLLEND SLLST X2REG2 ADD / 

ENGRID XI ISLOTEND RSLOT / 

DEFINE ISLOTEND ISLOTEND 1 ADD / 

STRUCTURE CONFORMAL SLOTRT -1 ISLOTST ISLOTEND SLRST SLREND 1 LEND / 

DO K 1 N 2 / 

DEFINE M1END LEND X3PER K MULTIPLY ADD / 

DEFINE MIST M1END X3REG1 MINUS X3REG1 MINUS / 

STRUCTURE CONFORMAL SLOTRT -I ISLOTST ISLOTEND SLLST SLLEND MIST Ml END / 

ENDDO / 

DO K 2 N 2 / 

DEFINE M1END LEND X3PER K MULTIPLY ADD / 

DEFINE MIST M1END X3REG1 MINUS X3REG1 MINUS / 

STRUCTURE CONFORMAL SLLT -1 ISLOTST ISLOTEND SLRST SLREND MIST M1END / 

ENDDO / 

COMMENT * FILL IN THE PORTION THE SLOT HOLLOWED OUT IN THE CONDUCTING WALL 


AT THE AXIAL END * / 

STRUCTURE CONFORMAL SLPATCH +1 1 X1TOTGRID 1 X2TOTGRID X3TOTGRID X3BLOCK / 

Q ************************ ******** ***★★★★****★***★***★★* **************** j 

Q a********************************************************************* j 

C GAP SECTION. . .GAP SECTON. . .GAP SECTION. . .GAP SECTION. .. GAP SECTION... / 

Q *★*■****★★★*★★★**★★■*★★★•***★★* + ★★***★★**★*★•*★★★***★**★*★**★•***★★*■*■★*★*★★ j 

COMMENT * FORM THE GAP * / 

DO K 0 N1 1 / 


DEFINE GAP 1ST X3REG1 X3REG2 ADD 1 ADD X3PER K MULTIPLY ADD / 

DEFINE GAP1END GAP1ST X3REG3 ADD / 

STRUCTURE CONFORMAL GAP1 -1 IFERRST IFERREN 1 X2TOTGRID GAP1ST GAP 1 END / 


ENDDO / 

Q **★★++*★★★**★+***+★****★**★**+******★*******★***★★*****+*+★*★★***★*+** j 

Q A********************************************************************* f 

C ft********************************************************************* / 

C DIAGNOSTIC SECTION DIAGNOSTIC SECTION .... DIAGNOSTIC SECTION / 

C ★**★★**★★★***★****★****★★**★***★★**★★★**★★***★★***★★*****★★★*★+★*★★★★* / 

GRAPHICS SYSTEM / 
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C DISPLAY MODEL 2 2/ 

C DISPLAY MODEL 2 12 / 

DIAGNOSE FREQUENCY 104/ 

FREQUENCY 4 150E+9 .001 50000 / 

RANGE 110. 00. 0 E311 11 1 X3TOTGRID / 

LINPRINT 1X3X1111 11 1 X3TOTGRID 001000000/ 

BALANCE 1 1 X1TOTGRID 1 X2TOTGRID 1 X3TOTGRID / 

OBSERVE 1 YES ENERGY 1 0.0 0.0 0.0 0.0 TOTAL 01/ 

c ********************************************************************** / 

Q **************** + ********************A****'ili***********A*************** f 

START / 

STOP / 
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